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Inelastic dark matter reconciles the DAMA anomaly with other null direct detection experiments 
and points to a non-minimal structure in the dark matter sector. In addition to the dominant 
inelastic interaction, dark matter scattering may have a subdominant elastic component. If these 
elastic interactions are suppressed at low momentum transfer, they will have similar nuclear recoil 
spectra to inelastic scattering events. While upcoming direct detection experiments will see strong 
signals from such models, they may not be able to unambiguously determine the presence of the 
subdominant elastic scattering from the recoil spectra alone. We show that directional detection 
experiments can separate elastic and inelastic scattering events and discover the underlying dynamics 
of dark matter models. 



I. INTRODUCTION 

The annual modulation anomaly from the DAMA 
experiment is an intriguing hint of the identity of dark 
matter [TJ [5]. Direct detection experiments such as 
DAMA measure the energy of nuclear recoils from in- 
cident dark matter particles. The signal experiences 
an annual modulation because of the variation of the 
Earth's relative motion with respect to the dark matter 
in the halo. While the DAMA experiment has measured 
an 8.2<t modulation with the correct phase, its results 
are in conflict with those from CDMS [H], XENON10 
0[7J, CRESST 015], ZEPLIN-II [ID], and ZEPLIN-III 
[llj . which measure the total, unmodulated scattering 
rate. 

All current direct detection experiments are opti- 
mized to look for elastic scattering, which has an expo- 
nentially falling recoil energy spectrum [12j . A distin- 
guishing feature of DAMA's measured modulation spec- 
trum is that it is suppressed at energies below ~ 25 
keVnr, where elastic events should dominate. DAMA's 
results, taken together with the null results from all other 
direct detection experiments, can be elegantly explained 
if the dark matter scatters inelastically off of nuclei to an 
adjacent state that is 0(100 keV) more massive than the 
ground state [TSHTB] . In the case of inelastic dark matter 
(iDM), a minimum velocity is needed to upscatter and 
the recoil spectrum is suppressed below this kinematic 
threshold. 

Inelastic dark matter has three features that lead to 
the consistency of all current experiments. First, iDM 
has a much larger dependence on the Earth's velocity 
than elastic dark matter and thus the annual modula- 
tion is often 25% or larger, in comparison to 2.5% for 
elastic scattering. This decreases DAMA's unmodulated 
cross section by a factor of 10 or more and therefore 
reduces the expected signal by the same factor in all ex- 
periments looking for unmodulated scattering. Second, 



iDM dominantly scatters off heavier nuclei and lighter 
nuclei may not have enough kinetic energy to excite the 
transition. As a result, CDMS' sensitivity is strongly 
suppressed due to inelastic kinematics. Lastly, inelastic 
scattering events have higher recoil energies than elas- 
tic events. Recent experiments such as XENON10 and 
ZEPLIN-III have shrunk their recoil energy window to 
eliminate higher energy scattering events, reducing the 
acceptance for iDM. However, XENON10 has recently 
reanalyzed their data over a larger signal window to be 
sensitive to iDM and their latest results are used in this 
article [7J. 

A dynamical alternative to the inelastic hypothesis 
was recently discussed in (TTJUH]; in this scenario, the in- 
teraction between the dark matter and Standard Model 
(SM) is elastic, but is suppressed by a form factor 

F dm (q 2 ) = c +c iq 2 +c 2 q i + ■■■ . (1) 

If the first term in this expansion vanishes, the elas- 
tic scattering rate is multiplied by additional factors of 
q 2 = 2m,NEn and goes to zero as Er — » 0. Form factor- 
suppressed scattering can arise from multipole, polar- 
ization, and charge-radius interactions between the dark 
sector and the Standard Model [19 . The recoil spec- 
trum for form factor elastic dark matter (FFeDM) is 
suppressed at low energies and more closely resembles 
the spectrum for inelastic events rather than standard 
clastic events, which peaks at low energies. However, 
it is more challenging to reconcile DAMA and the null 
experiments using only FFeDM because the form fac- 
tor depends on the product uinEr and any suppression 
for light nuclei can be compensated by looking at larger 
energies. 

In this paper, we consider a scenario where inelastic 
scattering is complemented by a form factor-suppressed 
elastic component. These scenarios arise, for example, in 
composite dark matter models |20j and yield new chal- 
lenges for direct detection phenomenology. As is illus- 



trated in Fig. [T] the modulation amplitude for inelas- 
tic scattering with subdominant charge-radius scatter- 
ing (dashed line) resembles that for inelastic scattering 
(solid line) and is markedly different from the regular 
elastic spectrum. Distinguishing the elastic and inelastic 
contributions is critical for understanding the underlying 
structure of the theory. This article addresses how di- 
rectional detection experiments, which can measure the 
direction of the nuclear recoil in addition to its energy 
[2T| [22] , can differentiate the contributions to the scat- 
tering rate. In the following section, we review the direct 
detection phenomenology. In Sect. |III[ we introduce di- 
rectional detection experiments and show how they can 
distinguish the dynamics in the dark sector. We con- 
clude in Sec. [TV] 



II. DIRECT DETECTION PHENOMENOLOGY 

Direct detection experiments measure the energy of 
nuclear recoils in a detector. For a detector consisting of 
nuclei with mass mjv, the differential scattering rate per 
unit detector mass is 



dR po I da 

= ( v 

dE R m dm m N \dE R 



(2) 



where pp = 0.3 GeV/cm 3 is the local dark matter den- 
sity. The rate depends on the dark matter mass m,j m , as 
well as the differential cross section. When both elastic 
and inelastic scattering are allowed, the differential cross 
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FIG. 1: Comparison of the modulation amplitude for two 
scattering scenarios: completely inelastic scattering (solid) 
and inelastic scattering with a subdominant elastic charge- 
radius component (dashed). The points show the modula- 
tion amplitude measured by DAM A and DAM A/LIBRA Q]. 
The electron equivalent energy (keVee) is the recoil energy 
rescaled by the quenching factor for iodine (qi = 0.085). 



section is parameterized by 
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where c e i,in are the dimensionless couplings for the elastic 
and inelastic interactions, A is the scale for new physics, 
and 



dap _ m N f7 N (f p Z + f n (A - Z)f 
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is the standard rate for elastic scattering off a nucleus 
with charge Z and atomic number A. Here, p is the 
reduced mass of the dark matter-nucleus system, on is 
the cross section for the dark matter-nucleus interaction 
at zero momentum transfer, and / Pjn are the couplings 
to the proton and neutron, respectively. Our results are 
normalized to / p = I and / n = 0, and we assume that 
A = I GeV. 

The scattering operators coherently couple the dark 
matter states to the nuclear charge, and the Helm form 
factor accounts for loss of coherence at large momentum 
transfer, g, 



\F B (E 



R)\ = 



3ji(|gko) 
\v\ro 



(5) 



where s = 1 fm, 



5s 2 and r = 1.2 A 1 / 3 fm 



The specific values of n m and n c \ are model- 
dependent. For instance, regular iDM has n in = 0, while 
regular elastic scattering corresponds to n e \ = 0. Form 
factor suppression in either the inelastic or elastic scat- 
tering interactions results in additional powers of E R . 
Composite inelastic dark matter models have n m = 1, 
though this does not make a substantive change to the 
spectrum as long as A 2 3> itinE r [2H ES]- In this pa- 
per, we focus on the scenario where inelastic scattering 
(n- m = 1) is supplemented by a form factor-suppressed 
clastic component (n c \ = 2), although all the conclusions 
are general and apply to n ul = 0. 

For a given dark matter velocity distribution func- 
tion f(v) defined in the galactic rest frame, 



da 
~d~E~R 



d 3 v f(v + v e ) v 
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~dE~n 



(G) 



The minimum velocity is set by the kinematics of the 
scattering process: 
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m N E R 



elastic 
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(7) 



where 5m is the dark matter mass splitting. The Earth's 
velocity in the galactic rest frame, v e , is defined as 



v e = v G +«e(i). 



(8) 
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In the coordinate system where x points towards the 
galactic center, y points in the direction of the galactic 
rotation, and z points towards the galactic north pole, 
the sum of the Sun's local Keplerian velocity and its 
peculiar velocity is |57] 



v Q (t) « (0, 220, 0) + (10, 5, 7) km/s. 



(9) 



The velocity of the Earth in the sun's rest frame is given 
by 



2n(t-t ) . . 2n(t-t ) 
ei cos h £2 sin 



yr 



yr 



, (io) 



where to is the spring equinox (w March 21), v® = 29. 
km/s is the orbital speed of the Earth [T2], and 



ei = (0.9931,0.1170,-0.01032) 
e 2 = (-0.0670,0.4927,-0.8676) 



(11) 



are the unit vectors in the direction of the Sun at the 
spring equinox and summer solstice, respectively [281I29] . 

Typically, the velocity distribution function f(v) is 
assumed to be isothermal and isotropic in the galactic 
frame. In the "Standard Halo Model" (SHM), the Kep- 
lerian velocity is constant throughout the galaxy and the 
velocity dispersion is assumed to be Gaussian. These as- 
sumptions are not consistent with recent N-body simu- 
lations of dark matter particles; simulations such as Via 
Lactea (30] show that the dark matter velocity dispersion 
is anisotropic and that the density falls off more steeply 
at larger radii than in the isothermal scenario. 

This article adopts the following ansatz for the dark 
matter halo velocity distribution 

f(v) = N(e~^ 2a - e-CWfo) 3 -) 6 (z; csc - »), (12) 

which reproduces the SHM in the limit a — ^ 1. This 
distribution function captures the most important qual- 
itative features of Via Lactea, but in an analytical form 
that is easy to compute with. The parameter a changes 
the shape of the profile near the escape velocity. Inelastic 
scattering events are particularly sensitive to the high- 
velocity tail because their minimum velocity is larger in 
comparison to the elastic case [3T] 131] ■ 

To find the regions of parameter space consistent 
with current direct detection experiments, we perform a 
global chi-squared analysis and marginalize over the six 
unknown parameters: vq, v csc , a, m^, 5m, and a p , the 
cross section per nucleon. The x 2 function is 

JVexp ( x pred — X ohs \ 
X 2 (mdm,Sm,a n ,v ,v csc ,a) = yW - " 



i=i 



(T, 



(13) 

where N CKp is the number of experiments included in 
the fit, Vf lcd is the predicted experimental result, X° hs 
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FIG. 2: The 95% contours in the md m — 8m parameter space 
corresponding to c c i/ci n = (blue), 0.4 (magenta), 0.6 (yel- 
low) and 0.8 (green) for the scenario with ni n = 1 and n e \ ~ 2. 
The four regions overlap one another. 



is the observed result and <7j is the known error. The 
velocity distribution parameters are constrained to be 
within 



200 km/s <v < 300 km/s 
500 km/s < v csc < 600 km/s 
0.8 < a < 1.25. 



(14) 



These values are motivated by rough observational con- 
straints [T2J [33] and the fact that (12 1 is a spherically 



symmetric form of the Via Lactea fits by [32] 

The results are fit to the first twelve bins (2-8 keVee) 
of the recoil spectrum measured by DAMA [T] . The mod- 
ulation amplitude from 8-14 keVee is combined into a sin- 
gle bin with amplitude -0.0002 ± .0014 cpd/kg/keVee. 
The contribution of this last bin to the total x 2 is typ- 
ically negligible. The region of parameter space that is 
consistent with DAMA is constrained by the null exper- 
iments, each of which has observed a number of events 
in its signal window. The results from CDMS, CRESST, 
ZEPLIN-II, and ZEPLIN-III are included in the fit, in 
addition to XENONlO's recently updated analysis with 
an expanded signal window [7] . We require that the the- 
ory not saturate the number of observed events for each 
experiment to 95% confidence. 

Figure [2] shows the allowed regions of m<jm — Sm 
parameter space for different ratios of clastic to inelastic 
scattering. The contours are defined as 



X 2 {mdm ,Sm,a n ,v , v csc , a) 



xLt + Ax^CL), (is) 



where A X 2 (95%) = 12.6 for six degrees of freedom. The 
minimal x 2 found was Xmin = 4.1 for c c \/c- m — and 
thus every prediction was forced to have a x 2 < 16.7 
corresponding to the 2a band. 

The best-fit point for the completely inelastic sce- 
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nano is 

{v ,v csc , a) = (279,586,0.80) (16) 
(m dm ,<5m,cr„) = (59 GeV, 121 keV, 5.9 x 10~ 39 cm 2 ), 

where <j n = a^{pb N / ji 2 ,) and /iN(n) is the reduced mass of 
the dark matter-nucleus (nucleon) system. As the ratio 
of elastic to inelastic scattering increases, the null ex- 
periments become more constraining and the minimum 
X 2 increases. The maximum allowed ratio of elastic to 
inelastic scattering is c c i/c- ln = 0.8. For this ratio, the 
best-fit point has \ 2 = 14-3 and corresponds to 



5.00 



(v ,v esc ,a) = (297,500,0.82) 
(m dm , 5m, a n ) = (74 GeV, 98 keV, 1.2 x 10~ 39 



(17) 



There is a high degree of degeneracy in the allowed val- 
ues of the halo parameters for each m^m ~ 5m point. 
The parameter points with the least tension have m,} m 
and 5m near the best-fit points of (16 1 and (17 1, but 



the halo parameters and cross sections can vary wildly. 
Some areas of the allowed regions are highly sensitive to 
the exact range of the halo parameters. For instance, 
the large mass region corresponds to exceedingly "slow" 
velocity distribution functions where vq < 225 km/s and 
a > 1.15. 

In the near future, definitive experiments such as 
XENON100 [34] and LUX [35] will confirm or refute the 
inelastic dark matter hypothesis. These two experiments 
are the upgrades to the XENON10 detector. The LUX 
detector will consist of 300 kg of liquid xenon with a 
planned fiducial volume of 100 kg, and it will be sen- 
sitive to recoil energies as large as ~ 300 keVnr [31)] . 
Figure [3] shows an estimated recoil spectrum at LUX af- 
ter 1000 kg-days and assuming 30% efficiency. The solid 
line denotes the spectrum for the best-fit point in the 
pure inelastic scenario, while the dashed line shows the 
spectrum for the best-fit point when c c i = 0.8cj n . LUX 
would see many events in its signal window; for models 
with a small elastic subcomponent, it could see as many 
as ~ 40 — 50 events in the winter (gray) and ~ 70 — 80 
events in the summer (black). 

The inelastic and elastic scattering channels have 
remarkably similar nuclear recoil spectra because of the 
threshold behavior of the dominant inelastic channel and 
the form factor suppression of low momentum transfer 
events. The shape of the spectra differ only at small 
recoil energies (< 10 keVnr). Because this is near the 
threshold energy of the experiment (~ 5 keVnr), it is 
difficult to unambiguously distinguish the two contribu- 
tions to the scattering rate with experiments such as 
LUX and XENON100. In the following section, we show 
that additional information about the nuclear recoil di- 
rection is a critical component in understanding the dy- 
namics of the dark sector. 
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FIG. 3: Estimated recoil spectrum at LUX with 1000 kg-days 
and 30% efficiency during the summer (black) and winter 
(gray). The solid line corresponds to the best-fit point for 
completely inelastic scattering; the dashed line corresponds 
to the best-fit point when c c i = 0.8ci n . 



III. DIRECTIONAL DETECTION 

Directional detection experiments take advantage of 
the daily modulation in the direction of the dark matter 
wind in the lab frame, which arises from the Earth's ro- 
tation around the galactic center [37J[3H]. I R particular, 
the direction of the dark matter changes every twelve 
hours as the Earth rotates about its axis. For the case 
of elastic scattering, the daily modulation amplitude can 
be nearly ~100%, compared to the ~2.5% change in the 
annual modulation amplitude. 

Measuring the strong angular dependence of the 
nuclear recoil can be an important tool for detecting 
dark matter [37] . Several directional detection experi- 
ments are currently running: DMTPC [32], New Age |40) . 
DRIFT 01], and MIMAC 05]. To get reasonable an- 
gular resolution, the track left by the recoiling nucleus 
must be sufficiently long (~ 1 mm). This means that 
the detector material must be a gas, because liquid and 
crystalline detectors have scattering lengths that are too 
long. All the current detectors are using CF4, except for 
DRIFT, which uses CS2. These detectors are specifically 
designed to look for elastic spin-dependent interactions, 
and so the use of atoms with high spin coupling, such as 
fluorine, is preferred. 

It was recently pointed out that directional detec- 
tion experiments can serve as important tests of inelastic 
dark matter [35] if much heavier atoms are used in the 
detectors. In this section, we show that such experiments 
are the key to distinguishing the contributions from elas- 
tic and inelastic scattering components in the dark mat- 
ter sector. To begin, we will briefly review how the rate 
equation derived in the previous section is generalized to 
include the angular dependence of the recoiling nucleus. 
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A more complete discussion of the theory can be found 
in [221113!. 

Consider the lab frame where the incoming dark 
matter particle has a velocity v' = vx and scatters off 
a nucleus at rest. The recoiling nucleus has velocity vr 
and makes an angle 9 with the x-axis. Energy and mo- 
mentum conservation yield an expression for the recoil 
velocity, 



Vr 



2fj,v 
m N 



cos 9, 



which can be written in the frame-invariant form 



V ■ V R 



0. 



V 



(18) 



(19) 



The differential directional scattering rate per unit de- 
tector mass is 



d 2 R 



Pa 



dEfjd cos 7 



(20) 

m dm m N ' •"- ' 
x / d 3 vf(v)d(v ■ v R -v e -v R - v min ). 



The integral expression for the differential rate is an ex- 
ample of a Radon transform, the properties of which are 
reviewed in |43j . The Radon transform for the modified 
halo distribution function we consider is 



f(w)= d 6 v f(v)5(v-v R -w) = 27r / dvvf(v), (21) 

J J w 

where w = v e ■ vr + v m - m = v e cos 7 + v m - ul . Note that 7 
is defined as the angle between the recoiling nucleus and 
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FIG. 4: cos 7 spectrum during the summer for the best-fit 
parameters when c c i = 0.8ci n , assuming a detector of CF3I. 
The spectrum is taken for E r — 50 keVnr. The dashed line is 
the total from form factor-suppressed and inelastic scatter- 
ing. The shaded circle marks the "cross-over" point where 
the elastic scattering starts to dominate over the inelastic 
scattering. 



the direction of the Earth's velocity v e 
integral for (12), we find 



Evaluating the 



f(w) = ttN{ {w 2 
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There are several important features of this expression. 
Firstly, the differential rate is peaked at cos 7 = — 1. this 
can also be seen explicitly from the rate equation, where 
the largest fraction of the parameter space is allowed 
when v e ■ vr — —1 in the delta function, which makes 
sense because the rate is maximized when the Earth is 
moving in the same direction as the dark matter wind. 

The second important feature arises from the theta 
function in f(w). In particular, there is a maximum 
value of cos 7 above which the rate is zero: 



COS7 n 



(22) 



This cutoff depends on the minimum velocity, which is 
larger for inelastic scattering than elastic scattering. As 
a result, the recoil spectrum for any inelastic scattering 
component will have a cut off below that of any clastic 
scattering component. This is illustrated in Fig.|4j which 
shows the recoil spectrum for the best-fit parameters in 
a c c i/ci n = 0.8 theory at a recoil energy of 50 keVnr, 
which is near the threshold of most current experiments. 
The rate for the inelastic interaction becomes negligible 
by cos 7 ~ —0.5. There is a tail at larger values of cos 7, 
which arises from the elastic scattering component. 

Whether the elastic and inelastic scattering compo- 
nents can be distinguished depends on whether there are 
enough events that fall along the tail of the spectrum. 
Fig.[5]shows the expected rate at the value of cos 7 where 
the transition from inelastic to elastic scattering occurs. 
The spread arises from varying over the six unknowns of 
both the particle physics and halo profile models. Even 
given the uncertainties of the halo profile distribution, 
one can obtain enough events on the cos 7 tail for dis- 
covery with approximately 5-10 kg years-worth of data. 
This is well within reach of current experiments; the 
dashed lines in Fig. [5] show the projected sensitivity of 
the DMTPC experiment, assuming a CF4 detector. If a 
heavier detector material is used, as would be needed to 
test for inelastic dark matter, then the sensitivity should 
be even larger. 



IV. DISCUSSION 

We have shown that directional detection experi- 
ments can distinguish mixed inelastic-elastic scattering 
scenarios that would otherwise be difficult to discover 
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FIG. 5: The detection rate corresponding to the value of cos 7 
at the "cross-over" point, where the elastic scattering starts 
to dominate over the inelastic scattering. The shaded regions 
correspond to c c i/ci n = 0.2 (magenta), 0.4 (yellow), and 0.8 
(green). The rate is calculated for summer at a recoil energy 
of E r = 50 keVnr, assuming a detector of CF3I. The dashed 
lines show the projected sensitivity for the current DMTPC 
CF 4 detector after 5 and 100 kg-yr [21]. 

unambiguously, even with next-generation experiments 
such as LUX and XENON100. In particular, directional 
detection experiments can detect elastic wide angle scat- 
ters that are kinematically forbidden in inelastic tran- 
sitions. To evade current null experiments, a relatively 
large exposure of several kg-yrs is needed, but this is well 
within reach of current detectors. However, these experi- 
ments are currently optimized to look for spin-dependent 
elastic scattering and need to use heavier nuclei, such as 
iodine or xenon, to have sensitivity to inelastic recoils. 

Distinguishing different scattering mechanisms is 
critical for understanding the underlying symmetry 
structure of dark matter theories. Many types of in- 
elastic models can give small elastic scattering contri- 
butions. For instance, nontrivial scattering mechanisms 
arise in models where the dark matter has a finite size, 



such as composite [2H P£MrJ] , atomic [S3] , mirror [37] , 
and quirky [48] dark matter. All these models have sev- 
eral scattering channels with a hierarchy of scales given 
by a series of higher dimensional operators. Both elastic 
and inelastic operators may be allowed, and approximate 
discrete symmetries may cause the elastic scattering rate 
to be subdominant to the inelastic rate [20 . 

Another class of models that leads to form factor- 
suppressed elastic scattering is characterized by a 
pseudo-Goldstone mediator between the dark sector and 
the SM [T7J. In this case, additional dimension six op- 
erators in the Lagrangian may no longer be negligible in 
comparison to the standard spin-independent and spin- 
dependent operators. These higher dimension operators 
lead to momentum suppression in the scattering rate. 
There has been much interest recently in models with 
light mediators [49H54] . given the results from PAMELA 
[55], Fermi [56], ATIC [57], and HESS [58]. 

If the dark sector does indeed have non-minimal 
structure, it may give rise to several types of scatter- 
ing mechanisms. Distinguishing these different scatter- 
ing events is of fundamental importance for understand- 
ing the symmetry structure and unraveling the dynamics 
of the dark sector. 
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